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ABSTRACT 

We have searched for optical identifications for 79 Chandra X-ray sources that lie within the half- 
mass radius of the nearby, core-collapsed globular cluster NGC 6397, using deep Hubble Space Tele- 
scope Advanced Camera for Surveys Wide Field Channel imaging in Ha, R, and B. Photometry of 
these images allows us to classify candidate counterparts based on color-magnitude diagram location. 
In addition to recovering nine previously detected cataclysmic variables (CVs), we have identified six 
additional faint CV candidates, a total of 42 active binaries (ABs), two millisecond pulsars (MSPs), 
one candidate active galactic nucleus, and one candidate interacting galaxy pair. Of the 79 sources, 
69 have a plausible optical counterpart. 

The 15 likely and possible CVs in NGC 6397 mostly fall into two groups: a brighter group of six 
for which the optical emission is dominated by contributions from the secondary and accretion disk, 
and a fainter group of seven for which the white dwarf dominates the optical emission. There are 
two possible transitional objects that lie between these groups. The faintest CVs likely lie near the 
minimum of the CV period distribution, where an accumulation is expected. The spatial distribution 
of the brighter CVs is much more centrally concentrated than those of the fainter CVs and the active 
binaries. This may represent the result of an evolutionary process in which CVs are produced by 
dynamical interactions, such as exchange reactions, near the cluster center and are scattered to larger 
orbital radii, over their lifetimes, as they age and become fainter. 

Subject headings: globular clusters: individual (NGC 6397) — X-rays: binaries — novae, cataclysmic 
variables 



1. INTRODUCTION 

Gravitational interactions involving "hard" binary 
stars (i.e. those with orbital velocities that exceed the 
local velocity dispersion) provide the energy source that 
drives globular cluster evolution from the time that the 



cluster core approaches collapse (Hut et al. 111992 Heggie 



& Hut 2003). These interactions have a profound effect 



on the long-term fates of both the cluster and the bina- 
ries. Among the hardest binaries in clusters are compact, 
accretional X-ray sources. Bright low-mass X-ray bina- 
ries (LMXBs), which contain a neutron star primary, are 
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strongly overabundant in dense clusters an d thus almost 



certa i nly have a dynamical origin there (see Pooley et al 



2003 Ivanova et al. 20081). There is also growing evi 



dence that a majority of the more numerous cataclysmic 
variables (CVs) in clusters may also have a dynamical 
origin, rather than re presenting a primordial population 
( |Pooley fc Hut|2006[ ). This is a complex issue, given the 
many possible channels for the formation of CVs in clus- 
ters, which include: the evolution of primordial binaries, 
mediated by hardening interactions; exchange interac- 
tions between singles and primordial binaries, possibly 
involving multiple exchanges; tidal captures which do 
not result in merger s; and three-body interactions (see 



Ivanova et al. 2006). Deep surveys of compact binary 
populations in globular clusters will help to clarify the 
relative importance of these CV formation mechanisms. 

NGC 6397 is the second clo sest globular cluster (d = 
2.3 kpc; |Strickler et al.|[2009[ ) and is by far the nearest 
cluster to have undergone core collapse. It thus provides 
a special laboratory for studying the formation and evo- 
lution of compact binary systems in the extreme envi- 
ronment of a collapsed core. NGC 6397 has been a key 
target for ROSAT, Hubble Space Telescope (HST), and 
Chandra X-ray Observa tory (Chandra) st udies of com- 



pact binaries in clusters. Cool et al. ( |1995[ ) first detected 



cataclysmic variables in JNGC 6397, as counterparts to 
faint ROSAT sources, using HST Wide-Field/Planetary 
Camera 1 Ha and i?-band ima ging to select Ha-excess 
objects. Grindlay et al.| ( |1995[ ) confirmed the power of 
the Ha— R photometric technique, providing HST Faint 
Object Spectrograph verifica tion of the first three CV 
identifications in NGC 6397. |Grindlay et~al| j2001| ) car- 
ried out the first Chandra observations of JNGC 6397, 
which detected a rich population of 25 s ources within 
the cluster half-mass radius (r h = 2.33'; |Harris| p96l) . 
WFPC2 imaging in Ha, R, V, and I has led to the lden- 
tification of many of these sources , including n i ne CVs 
and a number of active binaries (i Taylor l |2002[ [Grind 



lay||2006j ) . Since the study by Taylorj ( |20U2[ ) was carried 
out, we have obtained considerably deeper Chandra and 
HST imaging of NGC 6397. The results of the Chandra 
imaging, which detected 79 sources within the h alf-mass 
radius, are reported by Bogdanov et al. (20101. These 



X-ray sources, designated by a prefix U , are listed in Ta- 
ble [l] In this paper, we report a search for optical coun- 
terparts to these sources using a deep HST ACS/WFC 1 
dataset. We describe the data, the analysis method, and 
the results in the following sections. 

2. DATA 

The optical imaging used in this study is the HST 
GO-10257 dataset (PI: Anderson), which provides deep, 
highly dithered ACS/WFC imaging of the central region 
of NGC 6397 in F435W (B), F625W (R), and F658N 
(Ha) . One of the main motivations of this study was to 
search for wobbles in the positions of stars, indicative 
of massive unseen companions, such as white dwarfs, 
neutron stars, or black holes. To provide good time 
sampling, the center of NGC 6397 was imaged over 10 
single-orbit epochs, spaced at approximately one-month 
intervals between 2004 July and 2005 June. A by prod- 
uct of this observing strategy is that we imaged the 



cluster at many different orientations, allowing for an 
exquisite mitigation of detector-dependent photometric 
and PSF-related errors. In addition, the stacked images 
show round PSFs with no sign of diffraction spikes. The 
dataset consists of 5 short B (13s), 5 long B (340s), 5 
short R (10s), 5 long R (340s), and 40 Ha (390s) ex- 
posures. The short exposures were designed to provide 
accurate photometry for stars above the main-sequence 
turnoff (MSTO). Saturation only becomes an issue near 
the tip of the giant branch. With the large number of 
Ha frames, the PSF sampling is particularly good for 
this band. 

3. ANALYSIS METHOD 
3.1. Photometry 

The HST data were analyzed using software based on 
the program developed for th e ACS Globular Cluste r 
Treasury project, described in Anderson et al.| ( f2008a[ ). 
The routine first finds the astrometric and pnotomet- 
ric mapping from each exposure into a master reference 
frame. It then searches through the master reference 
frame one small patch at a time, identifying stars where 
a significant number of exposures indicate coincident de- 
tections. These stars are then measured two ways: (1) 
they are measured in each individual exposure where 
they could be found, and (2) they are also measured 
simultaneously in all the exposures (which provides a 
better measurement for the fainter stars). Stars were 
measured with spatially variable, library PSFs (see An- 
derson & King ACS ISR 2006-01 2 ) constructed from the 
GO-9444 data set (Pi-Cool) of u> Cen, using an aperture 
of whole pixels that was optimized for the brightness of 
the star and its particular surroundings. 

Since our best coverage was in Ha, we did our star- 
finding on those images, and then used the positions 
found to identify and measure the stars in the B and R 
images. The final photometry we report here comes from 
averaging the photometry from the individual stars, with 
sigma clipping applied to remove outlying values due to 
cosmic rays, defective pixels, etc. A total of 25004 stars 
were detected. For the B and R bands, the photome- 
try was performed independently for the short and long 
frames. Photometric c alibration to the VEG AMAG sys- 
tem was adopted from Strickler et al. (2009). 

In addition to the star finding and measuring, we also 
produced stacked images of the scene, in a procedure 
that is akin to the IRAF/STSDAS drizzle algorithm with 
pixfrac set to zero. The stacked images were over- 
sampled by a factor of two, in order to take advantage 
of the heavy dithering to increase the effective resolu- 
tion. The resulting supersampled stacked images have a 
12,000x12,000 format and cover an approximately circu- 
lar field of diameter 5' with a pixel scale of 0.025". Fig- 
ures [l] and [2] show a combination of the 40 Ha images, 
together with the error circles for the 79 Chandra sources. 
The error circles are plotted at their actual sizes in both 
figures. As can be seen in these figures, the combination 
of images at various orientations results in a very smooth 
symmetric PSF, with a distinct lack of diffraction spikes. 

Color-magnitude diagrams (CMDs) were constructed 
from the photometry using the R magnitudes, and the 



Advanced Camera for Surveys, Wide Field Channel 



2 http:/ /www. stsci.edu/hst/acs/documents/isrs/isr0601.pdf 
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Fig. 1.— Stacked HST ACS/WFC in the Ha filter of NGC 6397 
with Chandra source error circles. North is up and east is to the 
left. The source labels have been omitted for sources within 30" of 
the cluster center for clarity; these are shown in Fig. [2] The inner 
green circle represents the core radius of 5.5" and the outer green 
circle represents the half-mass radius of 2.33'. There are a total of 
79 sources detected within the half-mass radius. 
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Fig. 2. — The inner 1' X 1' region of the stacked Ha image with 
Chandra source error circles. The error circles are shown at their 
actual sizes. The green circle represents the core radius of 5.5". 



B — R and Ha — R color indices. The short-exposure B 
and R magnitudes were used for stars brighter than 
about 2.5 mag below the MSTO, since the long-exposure 
frames saturate above this point. Figures [3] and [3] show 
the resulting CMDs. The efficacy of the photometric 
procedure is indicated by the tightness of the fiducial se- 
quences. No proper-motion cleaning has been performed 
at this point in the a nalys is. Prop er-motion determi- 
nation is described in £3.2 In S 4.1 we will use proper 



motions in our source identification. Unfortunately the 
quality of the cluster-field separation depends too sen- 
sitively on brightness to provide a definitive member- 
ship determination for all stars. The (B — R, R) CMD 
reaches deepest for the bluest stars, since the faintest red 
main-sequence (MS) stars are below the detection limit 
in B. The upper part of the white dwarf (WD) cooling 
sequence is clearly detected in the (B — R, R) CMD, ex- 
tending to nearly 10 mag below the MSTO in R. There 
is also an indication of a second WD sequence above the 
primary sequence, which Strickler et al.| ( |2009 [) interpret 
as a He WD sequence. JNIote that in the (Ha— R, R) 
CMD, the WDs lie to the Ha-deficit side of the MS, re- 
flecting the strong Ha absorption lines in WD spectra rel- 
ative to those of faint MS stars. The blue stragglers and 
the blue horizontal branch stars similarly lie to the Ha- 
deficit sides of the subgiant branch and the giant branch, 
respectively. 

3.2. Astrometry 

The stacked HST ACS/WFC images for each filter 
were rectified to a common astrometric system. We cal- 
culated an absolute astrometric zeropoint for the stacked 
images relative to the ICRS using a two-step process. We 
first computed a plate solution for an ESO WFI frame 
(WFI.2002-06-21T01:45:08.153) using 357 primary astro- 
metric standards from the USNO UCAC2 catalog. We 
then selected 350 secondary astrometric standards from 
the WFI frame and used these to compute a plate solu- 
tion for our stacked ACS/WFC images. We determined a 
bores ight correction for the C handra source coordinates 
from Bogdanov et al. (2010) by computing the mean 
offsets between the HST and Chandra coordinates for 
the six brightest CVs and the millisecond pulsar (MSP) 
PSR J1740-5340. This resulted in a shift of the Chandra 
coordinates to the ICRS of -0.19" ± 0.02" in RA and 
0.18" ± 0.02" in Dec, where the quoted uncertainties are 
the standard deviations of the offsets for the seven ob- 
jects. As an external test of our boresight correction, we 
note that our optical position for PSR J 1740-5340 agrees 
with that of |Bassa fc Stappersj (I2004D to 0.01" in RA 
and 0.03" in Dec. They followed a similar astrometric 
procedure to what we used, with WFPC2 rather than 
ACS/WFC imaging. 

We searched for optical counterparts to the Chandra 
sources by overlaying the Chandra error circles on the 
HST stacked images, with the boresight cor rection ap- 
plied to the Chandra source positions from |Bogdanov| 
et al. (2010). The 9 5% confidenc e error circle sizes were 
computed following Hong et al. (2005). Since the un- 



certainty in the optical positions (^ 0.1") was small 
compared with the size of the X-ray error circle radii 
(~ 0.5"), we neglected the contribution of the former to 
the total positional uncertainty. 

A preliminary proper-motion analysis was carried out 
in order to do a first-pass cluster-field separation. The 
present dataset was used for the reference epoch, with the 
stellar positions averaged over the 10 orbits. Two differ- 
ent datasets were considered for the second epoch, GO- 
7335 (obtained with the WFPC2 in 1999) and GO-10775 
(obtained with the ACS/WFC in 2006). The former pro- 
vides a longer time baseline, while the later provides 
more depth and spatial coverage. To compute proper 
motions, we measured stars in each early-epoch expo- 
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Fig. 3. — Color-magnitude diagrams for stars within the half-mass radius of NGC 6397 and CV identifications. The new candidates have 
been selected based on their blue color and Ha excess. Note that in the right panel, the bright CVs lie to the Ho-excess side of the MS, 
while the faint CVs lie to the Ha-excess side of the WD clump, which itself lies to the Ho-deficit side of the MS. 



sure image, corrected these positions for distortion, then 
transformed these positions into the reference frame us- 
ing linear transformations based on only cluster mem- 
bers. By using only members to define the transfor- 
mations, we automatically get displacements (and hence 
motions) with respect to the cluster's systemic motion. 

The cluster-field separation that can be achieved for 
moderately bright stars is illustrated in Fig. [5j which 
shows the proper- motion distribution for stars with 18 < 
R < 19. Cluster and field regions are defined by a 
magnitude-dependent radius about the centroids of the 
apparent cluster and field clumps. Stars that do not 
reside in cither of these regions are considered to have 
an unknown proper-motion membership status. Since 
the proper-motion error increases with magnitude, this 
proper-motion separation becomes less effective with in- 
creasing magnitude. Eventually the cluster and field re- 
gions of the proper-motion plane show a strong over- 
lap. It was generally possible to obtain an approxi- 
mate discrimination between cluster and field stars to 
a limiting magnitude of R « 24. In a few cases 
where the 2006-second-epoch results were resulted in an 



unknown proper-motion membership status, the 1999- 
second-epoch results were definitive. 

4. RESULTS 

4.1. Chandra Source Identification 

Optical CMDs provide a powerful means of select- 
ing likely Chandra source counterparts and investigat- 
ing their properties. For each of the 79 Chandra sources 
within the half-mass radius, we checked the locations of 
all objects within the X-ray error circle in the CMDs 
shown in Figures [3] and HJ Objects that fell on the 
main sequence were considered to be unlikely counter- 
parts, given the relatively low X-ray to optical flux ratio, 
fx/fopt, of most MS stars, in contrast to the ranges for 
chromospherically active binaries and cataclysmic vari- 
ables. Table [I] summarizes the result of this counterpart 
search. Plausible identifications were obtained for nearly 
all sources and are listed in this table, along with pre- 
vious source designations, inferred object types, proper- 
motion information, and photometry. The only sources 
for which the error circle was empty are the quiescent 
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Fig. 4. — Color-magnitude diagrams for stars within the half-mass radius of NGC 6397 and AB/MSP identifications. The new candidates 
have been selected based on their red color and generally small Ha excess. The object U12 is a known MSP and U18 is a likely MSP. The 
object U70 is a blue straggler. Eight atypical ABs— U42, U62, U63, U65, U66, U86, U92, and U120— are considerably redder than typical 
and/or have larger than typical Ha excesses. 



low-m ass X-ray binary (qLMXB) U24 
200lj, U104, U106, Ulll, U113, U119, 



Grindlay et al. 
121, U122, and 

U123\ In the case of U24, there is a small "blip" near 



the center of the error circle in the stacked R image, but 
it likely represents the combination of Airy ring artifacts 
from two bright nearby stars. Several other sources, U5, 
U16, U77, and U84, lie outside of the R field and thus are 
difficult to classify, from these data, due to the lack of 
B — R and Ha — R color measurements. However, the U5 
and U77 counterparts, V30 and V36 respectively, have 
bee n shown to be yariabl e using ground-based ph otome- 



try ( [Kaluzny et~aLl|2006l ). 3 |Kaluzny et ah] ( |2006| ) classi 



fied V3U as an eclipsing binary that lies to the tar red of 
the main sequence (sug gestive of it being a non-member) . 
Rozyczka et al. (2010) have recently classified V36 as a 
single-line spectropsopic binary. We have classified both 
of these stars as uncertain active binaries (ABs). U5 
and U84 have very red B — Ha colors, suggesting that 
their B — R colors are likely very red as well, which sup- 

3 P refix V objects are variable sta rs fr om the ground-based s tud- 
ies of |Kaluzny fc Thompson] l |2003[ ) and |Kaluzny et al.|p006) . 



ports an AB identification. Table [T] also gives a sum- 
mary of the results of the proper-motion analysis, which 
indicates whether the object is consistent with being a 
cluster member, a field member, or neither, based on a 
comparison of its proper motion to the proper-motion 
distribution for stars of similar magnitude. 

4.2. Source Types 

Based on the location of the proposed counterparts in 
the CMDs (or, in a few cases, based on ground-based 
detection of variability) , we primarily assigned types of 
cataclysmic variable (CV) and chromospherically active 
binary (AB). CVs were defined as being significantly to 
the blue of the MS and having significant Ha excesses (ei- 
ther relative to the MS or to the white dwarf sequence). 
ABs were defined as lying within ~ 0.75 mag above the 
MS (and thus within ~ 0.2 mag to the red of the MS) 
and having small Ha excesses (< 0.1 mag) , based on the 
previo us analysis of ABs in NGC 6397 by |Taylor et al.| 
(2001 1. In three cases where the only object in the error 
circle was an apparent MS star, U41, U91, and U112, 
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Fig. 5. — Proper-motion components for those stars in the mag- 
nitude range 18 < R < 19 for which a proper-motion measurement 
was obtained. The 2006 dataset was used for the second epoch for 
this plot. The zeropoint corresponds, by design, to the systemic 
cluster motion. Stars that are clearly in the cluster distribution 
are plotted as green triangles and fall in a tight clump about the 
origin, stars in the field distribution are plotted as red squares, 
and stars with an unknown status are plotted as blue circles. The 
14 candidate identifications for stars in this magnitude range are 
indicated by blue numbers. Note that source U65 belongs to the 
unknown group and U63 does not have a measured proper motion, 
while the other sources in this magnitude range are consistent with 
cluster membership. A reasonable interpretation is that U65 is a 
relatively nearby active binary that is superposed on the cluster; 
see Fig. H] for the photometric status of this star. 



we note its presence in Table [T] Similarly, an apparent 
MSTO star is the only object present in the error circle 
for source U118. We note that an AB with a low mass 
ratio and weak lines could look like a MS star in both 
CMDs. 

Two objects, U28 and U108, were classified as back- 
ground galaxies based on the extended appearance of 
their images. U28 resembles an edge-on spiral, while 
U108 has a more complex structure suggestive of inter- 
acting gal axie s. The high X-ray to optical flux ratio of 
U28 (see £4.3) suggests that it is an active galactic nu- 
cleus ( AGJN ) . While the flux ratio of U108 was not deter- 
mined, due to the complex nature of the optical image, 
it is clearly much lower than that of U28. One possibility 
is that we are detecting X-ray emission from U108 that 
is produced by a galaxy collision. 

Figure [3] shows the location of the CV candidates in the 
CMDs. There is a suggestion of an evolutionary sequence 
from young, bright CVs to old, faint ones. We return 
to this point in ^JEJ The six brightest CVs mostly lie 
about 0.2 - 0.8 mag to the blue of the MS. The optical 
emission of the systems appears to be dominated by the 
secondary in the R band, with a larger contribution from 
the disk in the B band. The fairly high i?-band flux 
indicates that the secondaries are relatively ma ssive, ~ 
0.5 — .7 MfT >, as inferred from the isochrones of |Baraffe 
et al.| (|1997[). All of the bright CVs have substantial 



Ha excesses relative to the MS; these excesses generally 
increase with magnitude. 

The CV candidates 1-8 are numbered accord ing to the 
schem e used in previous papers, as listed by |Grindlay 



et al. (2001). A prefix CV is used to designate each 



CV candidate. |Grindlay et al.| ( |2001[ ) identify U28 as a 
CV; however we find that this Chandra source appears 
to be a galaxy. Here we use the identification by |Taylor| 
p002| of U60 as CV 9. CVs 10-15 a re newly identified 
in the present study; Grindlay (2006) assigned numbers 
for CVs 10 and 11. We have numbered CVs 12-15 based 
on optical luminosity within this group. The translation 
between Chandra source numbers and CV numbers is 
given in Table [l] 

There is a two magnitude gap between the bright CVs 
and the fainter CVs. Just below this gap are two possible 
transitional objects, U7 and U61, which lie between the 
white dwarf sequence and the main sequence. We note 
that the proper-motion data produce an unknown mem- 
bership status for U61, while U7 registers as consistent 
with membership. This leaves open the possibility that 
U61 is a background AGN with an emission line that is 
shifted into the Ha window. However, the image of the 
U61 counterpart does not show obvious evidence of ex- 
tension. The fainter CV candidates, below U7 and U61, 
lie in the vicinity of the WD cooling sequence. The op- 
tical fluxes for the faint CV candidates are clearly domi- 
nated by the contribution of the WD. The Ha— R indices 
for these faint CVs nearly all lie at least 0.5 mag to the 
Ha excess side of the WD clump, which is centered at 
about Ha— R = 3.25. This suggests that the faint CVs 
have a strong Ha-emission core (due to an accretion disk) 
within the broad absorption lines of the WD continuum. 
This inferred spectroscopic property of the faint CVs ap- 
pears to be generally consistent with th ose of the WZ 



Sge c lass of evolved cataclysmic variable ( Schwarz et al. 
2004|. Another possibility is that some 01 these objects 



are AM CVn type double-degenerate systems, in which 
a low-mass He WD donor feeds a much more massive 
carbon-oxygen WD. Such systems show no evidence of 
hydrogen lines in their spectra and instead have He I 
and sometimes He II lines, which are typically in emis- 
sion ( |Anderson et al.|2005||2008b[ ). AM CVn stars could 
have evolved into contact from the known He WDs in 
NGC 6397, which have bee n shown to have (det ached) 
heavy CO WD companions ( |Strickler et al.||2009| ). Such 
an object should lie close to the main sequence in the 
(Ra-R,R) CMD. Possible candidate AM CVn stars 
among the CV candidates reported here are U25, U31, 
U60, U80, and U83. 

We note that most of the WD-like CV candidates clus- 
ter near a magnitude of R ~ 23.5, which corresponds to 
usi ng the distance modulus determined 
20091 



a M 



11.2, 



by Strickler et al 



We note that Gansicke et al. 
( |2009p have reported a sharp peak ("spike") near the 
period minimum (P ss 80 — 86 min) in the period distri- 
bution of CVs selected from the Sloan Digital Sky Survey 
(SDSS). The average absolute magnitude for the CVs in 
the spike is (M g ) = 11.6 ± 0.7. Allowing for the differ- 
ence in photometric band, this is reasonably close to the 
characteristic magnitude of the faint CVs in our study. 
This suggests that these faint CVs may well belong to 
the spike population, i.e. have periods in the vicinity of 
the period minimum. This could be tested by deter- 
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mining orbital periods for the faint CVs fr o m su itable 
time-resolved photometry. Gansicke et al. ( |2009| find 
that of 33 CVs with orbital periods less than 8b min, 
20 have spectra that are WD-dominated. In our sample, 
it appears that all seven faint CVs are WD-dominatcd. 
Given the small sample size, it is not clear whether this 
represents a significant difference between faint field and 
faint cluster CVs. 

Figure [4] shows the location of the objects identified 
as AB stars in the CMDs. We note that these objects 
sometimes lie near the edge of the MS in either the right 
or left panel, but deviate by a larger amount in the other 
panel. Of the 42 AB counterparts listed in Table[T] 26 are 
newly identified in our study. Overall, the ABs appear 
to form a relatively homogeneous binary sequence along- 
side the MS. presumably mostly differing in mass. [Taylor 
et al. (2001) have argued that these stars are likely BY 
Draconis stars and represent a hard binary population in 
NGC 639 7. In a large study o f photometric variability 
in 47 Tuc, Albrow et al. (2001) found a substantial pop- 
ulation of BY Draconis binaries, as well as W UMa and 
other contact binaries, and eclipsing binaries. 

In addition to the likely BY Dra population, there are 
a number of counterparts that show a significantly dif- 
ferent distribution in the CMD. These include the "red 
straggler" counterpart to U12 (the millisecond pulsar) 
and the similar counterpart to U18, the blue-straggler 
counterpart to U70, a group of five stars with very red 
colors, and a group of five stars with strong Ha excesses. 
There are two stars in common between these latter two 
groups. 

The counterparts to U12 and U18 both have proper 
motions consistent with cluster membership. Given their 
similar optical and X-ray properties, it a ppears likely 
that UI8 is a MSP flBogdanov et al.||2010[ ). The coun- 
terpart to U70, which registers as a proper-motion clus- 
ter member, has been identifi ed as variable star V20 by 



Kaluzny fc Thompsonj ( |2003[ ). |Rozyczka et al.| ( |2010 | 
find that it is a double-line spectroscopic binary with a 
period of 0.86 d and a mass ratio of q ~ 0.2. The He WD 
candidate PC-5 ( Taylor et al.|200T ), which has been pre- 
viously suggested as a possible counterpart to U70, lies 
well outside of the current error circle. 

The counterparts to U42, U63, U65, U86, and U92 dif- 
fer strikingly from the typical ABs, lying up to 1 mag to 
the red of the MS, and thus well above it. Of these, the 
counterparts to U42 and U92 hav e been previously iden - 
tified as optically variab le st ars ( |Kaluzny et al. 2006), 
as confirmed here (see ^4.4 and Fig. [7]); the counter- 
parts to U63, U65, and U86 are new identifications. The 
U9 2 counterpart (V7) is c l assifie d as a W UMa system 
by |Kaluzny fc Thompso"n| ( |2003| , while the U42 coun- 
terpart ( V 26 ) i s clas sified as an irregular variable by 
Kaluzny et al.| ( |2006[ ) . U92 is consistent with the clus- 
ter proper-motion distribution. The proper-motion in- 
formation available for U42, U65, and U86 produce an 
unknown membership status. These stars have proper 
motions that are well outside of the cluster distribution 
but they are also inconsistent with the field distribution. 
There is no proper-motion information for U63. One 
reasonable possibility is that these stars are foreground 
active binaries superp osed on the cluster, as discussed by 
Kaluzny et al] ((2006} in the case of U5 (V30). In this in 



dwarfs with distance moduli that put them well above 
the cluster main sequence in the (B — R, R) CMD. 

Figure [4] also shows several AB candidates that have 
larger than typical Ha excesses; these include the coun- 
terparts to U62, U65, U66, U86, and U120. There is 
no proper-motion information available for U62. The 
proper-motion information for U65, U86, and U120 pro- 
duces an unknown membership status. U66 is consis- 
tent with the field proper-motion distribution. Thus, 
the proper-motion information does not offer much assis- 
tance in determining whether these stars are likely cluster 
members. 



terpretation, the foreground binaries contain K/M type 



Rozyczka et al. (2010) have drawn attention to two of 
the variable stars that we have identified as AB coun- 
terparts, V17 and V36, which correspond to the sources 
U76 and U77. These similar objects both lie close to the 
main-sequence tu rnoff in the CMD presented by |Rozy-| 
czka et al. (20101. Both objects are single-line spectro- 
scopic binaries, m which the unseen primary is inferred to 
have a mass in excess of 1 M . They suggest that these 
two objects are members of a class of dormant degen- 
erate binaries which may have much larger peak X-ray 
luminosities. 

4.3. X-ray to Optical Flux Ratios 

As an additional aid to the classification of sources, we 
have examined the X-ray to opti cal flux ratio, fx /foot , 
where we take fx (0.5—6 keV) from Bogdanov et al. ( 2010 1 
and set / opt = f R = 1.52 x 10 -o.4K-b The latt g r con _ 
version factor is computed from the HST flux calibra- 
tion constants. The X-ray count to flux transformation 
assumes a power-law spectrum with a photon index of 
T = 2.5 ( [Bogdanov et al ||20 10|) The resulting flux ratio 
is plotted versus fx in Fig. [B| The ratio is observed to be 
higher for accretional sources, such as CVs, LMXBs, and 
AGN s, than for chromospherically active binaries (see 
e.g. |Bassa et al.||2004[ who studied the Chandra source 
distribution in M4). Consistent with this, Fig. [6] shows 
that the CV and AB candidates form two distinct groups 
with no overlap in flux ratio. The median flux ratio is 
about 500 times larger for the CVs than for the ABs. 
Within each group there is a broad range of the flux ra- 
tio, about a factor of 40 for the CVs and a factor of 100 
for the ABs. There is no apparent dependence of the 
flux ratio on fx over a range of more than three decades 
in fx for the CVs and nearly two decades for the ABs. 
The apparent AGN, U28, has the highest flux ratio of 
all objects. The known millisecond pulsar, U12, has a 
higher flux ratio than most of the ABs and is about 40 x 
brighter in the X-ray than a typical AB. The object U18, 
which has a similar optical counterpart to that of U12, 
has even higher values of fx and fx / /opt ■ This supports 
the interpretation of U18 as a MSP not yet detected at 
radio wavelengths (Bogd anov et al.||2010[ ). 

Interestingly, the X-ray flux of CV5 (U22) is consider- 
ably lower in the 2007 observations than in the previous 
2002 observations; the flux differs more than an order of 
magnitude. This gives CV5 the lowest fx/fopt value of 
all the CVs when the 2007 X-ray flux is used together 
with the 2004-05 optical flux. Since the X-ray and op- 
tical fluxes were not measured contemporaneously, it is 
possible that CV5 typically has a flux ratio that is more 
consistent with the other bright CVs. We note that the 
X-ray flux of a bright CV typically drops during out- 
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Fig. 6. — X-ray to optical ij-band flux ratio vs. X-ray flux (0.5- 
6 keV) for CVs (blue triangles), ABs (red squares), MSPs (green 
diamonds), unidentified objects and main sequence stars (black 
squares), and an AGN (inverted magenta triangle). The ij-band 
flux is taken as a measure of the optical flux. The upper axis gives 
the equivalent X-ray luminosity assuming that all objects are at 
the distance of the cluster. Note that the ABs lie entirely below 
the CVs. 



burst (e.g. Wheatley et al. 2003), suggesting that CV5 
may have been in outburst in 2007. The objects U62, 
U66, and U86, which lie among the ABs with the four 
highest flux ratios, all have significantly larger Ha ex- 
cesses than the typical AB. It is not surprising that two 
measures of chromospheric activity, X-ray flux and Ha 
flux, are correlated for these stars. 

4.4. Variability 

Since our dataset provides a 40-exposure time sequence 
of Ha exposures, with four exposures per orbit, it was 
possible to investigate optical variability. The time se- 
quence samples time scales shorter than about one hour 
(the visibility period per HST orbit) and also time scales 
from one month to one year. The analysis of variability 
is complicated by the presence of outliers in the time 
sequences, some of which represent photometric anoma- 
lies. We investigated several measures of variability and 
adopted the RMS deviation about the mean of the Ha 
magnitude measurements, computed using an iterative 
3(7 clip to reduce the impact of outliers. We plot <r(Ha) 
versus mean Ha magnitude in Figure [7J This procedure 
results in a measure of variability that is most sensitive 
to orbital variability of binary systems, rather than large- 
amplitude fluctuations of CVs, given the outlier filtering. 
Thus, we also investigated the total range of the Ha mag- 
nitudes for each object and constructed light curves for 
each of the CV candidates in order to test for outburst 
behavior. 

As can be seen in Fig. [7J most stars fall on a well- 
defined sequence of increasing cr(Ha) with increasing 
magnitude. Stars that lie significantly above this se- 
quence generally show evidence for variability, although 
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Fig. 7. — Ha variability versus Ha magnitude. The ordinate is 
the dispersion of the up to 40 Ho measurements for each star. A 3<r 
clip has been used to filter out photometric problems. The set of 
stars that registers as variable, i.e. lies above the main distribution, 
is mostly independent of whether or not sigma clipping is applied. 
Sources with at least nine measured magnitudes are plotted. 



some of this scatter is likely due to photometric issues. 
We note that the mean trend of the variability index 
cr (Ha) with mean Ha magnitude differs for stars brighter 
than Ha ss 18, for which cr(Ha) increases from a min- 
imum of about 0.003 mag to about 0.04 mag for the 
brightest stars plotted. This behavior for bright stars 
is due to saturation issues that increase in size with de- 
creasing magnitude. 

We have plotted the positions of the previous and new 
Chandra source counterparts in Fig. [7j Many of these 
counterparts show significant variability as measured by 
the a of the Ha time series. This group includes all of 
the bright CVs 1-6: sources U10, U17, U19, U21, U22, 
and U23. Of the newly identified faint CV counterparts, 
only U61 registers as clearly variable by this measure and 
U25 shows a possible slight variability. While the other 
faint CV candidates have cx(Ha) values that are compa- 
rable in size to those of the bright CVs, these values are 
not significantly larger than those of other stars of sim- 
ilar magnitude. Thus, the failure to detect statistically 
significant variability for the other faint CV counterparts 
may well be due to the decreasing sensitivity of the vari- 
ability test with increasing magnitude. 

The total magnitude range for each of the six bright 
CVs was about 0.6-1.3 mag, which is less tha n the 2- 
5 ma g amplitude for a dwarf nova outburst ( jWarner 
1995| | . Source U21 (CV4) had the largest amplitude (1.3 
mag) of these bright CVs, with a light curve that a high 
phase of about 50 days duration. Source U61 (CV12) 
displayed a total range of 2.3 mag. The faintest mag- 



nitude measurements, which contribute significantly to 
this large range, have a high degree of uncertainty. The 
light curve shows a clearly differentiated low state in the 
first half of the year of observation and a high state dur- 
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ing the second half. The mean offset between these two 
states is about 1 mag and thus also falls short of a dwarf 
nova outburst. 

A number of the candidate AB counterparts also show 
evidence of variability, including six newly identified 
ABs: sources U79, U97, and U114, which show strong 
variability at the level of the six bright CVs; sources U88 
and U101, which show moderate variability; and source 
U103, which shows slight evidence of variability. Five 
previously identified AB counterparts, U14, U18, U42, 
U70, and U92, show clear variability, and U76 and U81 
show weaker evidence of variability. We note that U18 
has a similar level of variability to the MSP counterpart 
U12, again suggesting that these are similar objects. 

We carried out a period-folding light curve analysis for 
the bright CVs 1-6 and the MSP U12. Given the very 
unequal spacing of the data, it is difficult to determine 
a period ab initio, due to the high degree of aliasing. 
However, we were able to confir m the previously deter- 
mined orbi tal periods ob tained by |Kaluzny fc Thompson 
(|2003b and |Taylor| 02002] for C Vs 1 and 6, as well as that 
of U12 ( |D'Amico et al.| 2001). We also see evidence for 
additional long-term variability for CVs 1 and 6, at the 
level of a few tenths of a magnitude. 

4.5. Spatial Distribution 

We determined the cluster center by iterative centroid- 
ing in a 20" radius aperture using a sample of MS stars 
with magnitudes in the range 16 < R < 22. The result- 
ing center of a = 17 h 40 m 42.17 s , S = -53° 40' 28.6" 
lies within about 1" of CV2 (U19). Experimentation 
with the centroiding aperture size and the stellar sample 
definition indicates that the center position is uncertain 
by about 1". We then determined the cumulative radial 
distributions of a number of stellar groups out to a ra- 
dius of 100", which is the approxima te size of the centra l 
surface density cusp of NGC 6397 ( |Lugger et al. 1995). 
Figure M shows the distribution functions tor the stellar 
groups listed in the figure legend. As has been previ- 
ously noted, the bright CVs in NGC 6397 show a strong 
central concentration relative to other stellar groups, in- 
cluding the fainter CVs. With the exception of CV6, the 
five other bright CVs lie within 11" (2r c ) of the cluster 
center. In comparison, with the exception of U25, the 
faint CVs lie at least 20" (4r c ) from the cluster center. 
We carried out Kolmogorov-Smirnov intercomparisons of 
the stellar samples defined in Fig. M The comparison of 
each group to a MSTO sample, defined by 16 < R < 17, 
is given in Table [l] where the probability, p, of the two 
samples being drawn from the same parent distribution 
is listed. The bright CV and blue straggler (BS) sam- 
ples both differ very significantly from the MSTO sample 
(p < 0.1%), while the overall CV and AB samples differ 
from the MSTO sample at the 4% level and 1% level, re- 
spectively. The bright and faint CV samples differ from 
each other at a significance level of 3% and the bright 
CV and AB samples differ at a level of 2%. 

In order to further investigate the spatial distribution 
of the Chandra sources in NGC 6397 and the implications 
for object masses, we carried out maximum likelihood 
fits of what has been termed a "generalized King model" 
to the surface density distributions. This model can be 
described more generally as a "cored power law," in that 
it has a smooth transition from an inner core region to 
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Fig. 8. — Cumulative radial distributions for selected stellar 
groups. For each group the actual distribution (stair-step curve) 
and maximum likelihood fit (smooth curve) are shown. Note that 
the bright CVs show the highest degree of central concentration, 
followed by the blue stragglers. The distribution of faint CVs is 
statistically consistent with the MSTO stars. Fitting information 
for these stellar groups is given in Table ^ 



an outer power law. It takes the form, 



S(r) = S 



a/2 



(1) 



with the core radius r c related to the scale parameter r 
by, 

/ , \ !/ 2 

r c = [2- 2 ' a - l) r . (2) 

An "analytic King model," which provides a good fit 
to the surface density profile of the turnoff-mass stel- 



lar group fo r a normal-core cluster (e.g. 47 Tuc; Heinke 
|et al.|[2005| , corresponds to a = — 2 and r c = r . Since 
1NGC 6397 has a central-cusp structure, due to its core- 
collapsed status, we expect that the turnoff-mass group 
will h ave a flatter power-law slope than a = —2 (Colm 
1985 1. We performed a maximum- likelihood fit of Fqn.ll 
to a turnoff-mass group, defined byl6<i?<17 ana 
radial offset within 100" of the cluster center, to deter- 
mine the best-fit values of a and r e . We used nonlin- 
ear optimization to maximize the likelihood and boot- 
strap resampling to estimate the paramet er uncertain- 
ties, as discussed by Grindlay et al. (2002). The result- 
ing best-fit parameter values are a = —0.93 ± 0.05 and 
r c = 5.5" ± 3.4". 

By fitting Eqn.[l]to the surface density distributions of 
individual groups such as the CVs and ABs, the charac- 
teristic m asses of these object s can be estimated , as de - 



scribed by Heinke et al. ( 2005 ) and Lugger et al. ( 2007 ) 



In order to do tins, we adopt the approximation here 
that stellar groups with masses exceeding the turnoff 
mass are in approximate thermal equilibrium with the 
turnoff mass group and that the velocity dispersion of 
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TABLE 2 

Cored Power-Law Model Fits to 100 Arcsec 



Sample 


N a 


<? 


rc (") 


ct 


m (M©) 




K-S prob c 


MSTO 


1111 


1.0 


5.5 ±3.4 


— 0.93 ± 0.05 


0.80 ± 0.05 






CV 


14 


1.42 ±0.16 


3.3 ±2.1 


— 1.74 ± 0.32 


1.14 ±0.14 


2.3 


4.0% 


bright CV 


6 


1.83 ±0.26 


2.5 ± 1.6 


-2.53 ±0.51 


1.46 ±0.22 


2.9 


0.07% 


faint CV 


8 


1.04 ±0.24 


5.1 ± 3.8 


-1.01 ±0.47 


0.83 ±0.20 


0.2 


89% 


AB 


36 


1.32 ±0.08 


3.6 ±2.2 


-1.54 ±0.17 


1.06 ±0.08 


2.6 


1.1% 


BS 


23 


1.52 ±0.10 


3.0 ±1.9 


-1.93 ±0.21 


1.22 ±0.10 


3.7 


0.03% 



a Size of sample within 100" of cluster center 
b Significance of mass excess above MSTO mass 
c K-S probability of consistency with MSTO group 



each mass group is approximately constant with radius 
within the central cusp region. This approach is moti- 
vated by Fokker-Planck simulations of post-collapse core 
oscillations which indicate that a cluster core spends a 
majority of the time near a maximally expanded state 
in which the more massive stellar group s are close to 
thermal equilibrium (Murphy et al. 1990). More recent 



Monte-Carlo and N-body simulations suggest that the 
presence of primordial binaries may delay the gravother- 
mal oscillation phase to well beyond a Hubble time and 
that even clusters with an observationally core-collapsed 
designation suc h as NGC 639 7 may still be in the binary 
burning phase (Fregeau 1[2008| ). In any case, the assump- 
tion of thermal equilibrium among the most massive com- 
ponents in the central region of an evolved cluster is a 
useful starting point. 

In the thermal equilibrium approximation, the surface 
density profile for a mass group with mass m is given by 
Eqn. [I] with a slope parameter a related to the turnoff- 
mass slope a to by 



q (a to - 1) + 1 



(3) 



where q — m/m to . Thus, fitting this model involves just 
one free parameter, the mass ratio q, as So is determined 
by normalization and tq and ato are determined by the 
fit to the turnoff-mass group. 

Table [T] gives the results of maximum-likelihood fits 
of Eqn. [lj to the turnoff mass stars, CVs, ABs, and a 
comparison group of blue stragglers (BSs). The cumula- 
tive radial distributions for these groups and the model 
fits are shown in Fig. [HJ As can be seen from the ta- 
ble, the q values for all of the groups exceed unity, indi- 
cating that the characteristic masses exceed the turnoff 
mass. For all cases other than the faint CVs, the ex- 
cesses are significant at the 2a or higher level. The groups 
with the strongest evidence for a mass significantly above 
the turnoff mass are the bright CVs (2.9a) and the BSs 
(3.7a). We note that the results of this analysis agree 
with the K-S comparison results given in the last column 
of the table. The inferred mass range for the bright CVs 
(1.5 ± 0.2Mq) is consistent with a system made up of 
a heavy white dwarf (e.g. Mwd ~ O.9M0) with a main 
sequence secondary of mass ~ 0.6M Q . The inferred mass 
range for the faint CVs (0.8 ± 0.2M Q ) is consistent with 
a somewhat less massive white dwarf (e.g ~ 0.7 M Q ) 
and a secondary mass that has been whittled down to 
~ 0.1M Q . In such a system, the optical flux would be 
dominated by the white dwarf, as observed here. We pre- 
sume that the white dwarf in this case has been heated by 



the accretion process, since noninteracting heavy white 
dwarfs are too faint to detect due to the degree of cooling 
that has occurred over their long lifetimes. As discussed 
in §4.2| an alternative interpretation is that some of the 
CVs are AM CVn systems, in which the secondary is a 
He WD. In this case, the mass of the secondary would 
likely be about 0.2-0.3 M ( |Strickler et al.|[2009[ ). It is 
interesting to compare the wh ite dwart masses inferred 
here with those for field CVs. |Southworth et al.| (12009 1 
have recently determined a mass of 0.78±0.12 Mq for the 



WD in a field C V that lies above the period gap. |Little- 
fair et al. ( 2008 ) have determined the masses of the WDs 
in seven CVs that lie below period gap, with a resulting 
mean mass and standard deviation of 0.87 ± 0.07 M . 
Thus, given the uncertainty ranges in the CV masses in- 
ferred here, our results for likely WD masses are generally 
consistent with these measured WD masses. 

5. SUMMARY 

We have shown that with moderately deep B, R, and 
Ha ACS/WFC imaging, it is possible to identify 

and classify nearly all of the 79 currently detected Chan- 
dra sources within the half-mass radius of NGC 6397. 
Our classifications are based on CMD location and are 
tested for consistency with expected ranges of the X-ray 
to optical flux ratio. In only a few cases was it not possi- 
ble to classify the source, generally due to a lack of cover- 
age in one filter. The vast majority of the sources in the 
cluster can be classified as either cataclysmic variables (a 
total of 15) or active binaries (a total of 42). The X-ray 
to optical flux ratio for the CVs (0.2 - 6) substantially 
exceeds that for most of the ABs (0.0003-0.2). There is 
one previously known MSP and one potential MSP; the 
latter classification is based on the similar photometric 
properties of the two objects. 

The proper-motion test for cluster membership pro- 
duces an unknown status for many of the Chandra source 
counterparts. This is largely a consequence of the short 
time baseline of 1.5 years between the two epochs of 
ACS/WFC data used for determining proper motions. 
While the 1999 WFPC2 observations give a longer base- 
line, fewer stars are detected by the WFPC2. Thus, the 
proper-motion part of this study would benefit substan- 
tially from new HST observations of NGC 6397 with 
either the ACS/WFC or the WFC3. While it appears 
likely that the vast majority of the proposed counter- 
parts are cluster members, proper-motion confirmation 
would be valuable. It would be particularly interesting 
to know how many of the atypical ABs are cluster mem- 
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bers. While it is likely that most are foreground stars, 
one of the very red ones has appears to be a member 
on the basis of our preliminary proper-motion measure- 
ments. 

The CV distribution in NGC 6397 has a bimodal char- 
acter. We detect distinct bright and faint populations 
with only two possible transition objects. The optical 
emission of the six brightest CVs appears to be domi- 
nated by a relatively massive secondary, while that of 
the faint CVs appears to be dominated by the WD, with 
very little contribution from a very low-mass secondary. 
Alternatively, some of the faint CVs may be AM CVn 
systems, with a He WD secondary. The distribution of 
CVs in optical luminosity and color is consistent with ex- 
pected CV evolution due to the effects of mass transfer 
and angular momentum loss. As CVs age and the sec- 
ondary loses mass to the primary, the orbit tightens due 
to magnetic braking and ultimately gravitational wave 
radiation. In this picture, the bright CVs 1-6 represent 
young, recently formed systems, while the faint CVs rep- 
resent old, highly evolved systems. While CVs can be 
born at any orbital period, depending on the nature of 
the secondary star, there is a bias for CVs formed in 
exchange interactions in collapsed cluster cores to have 
more massive secondaries and thus longer orbital peri- 
ods than newly formed field CVs. This is both the result 
of mass segregation in the cluster, which may invert the 
mass function in the cluster core, and of the tendency 
for more massive stars to displace less massive ones in 
exchange interactions. Thus, while it is possible for a 
faint CV to be a young system, we expect that the bulk 
of the faint CVs evolved from an earlier bright state. 

Given the likely ongoing formation of CVs by dynam- 



ical processes in the central region of the cluster, it is 
plausible that we are observing a roughly equilibrium 
CV population in NGC 6397. However, there remains a 
possibility that the current population of six bright CVs 
represents the result of a recent core collapse event that 
resulted in an enhanced production of CVs. In this case, 
the bright CVs would be overpopulous relative to the 
faint CVs. 

The six brightest CVs in NGC 6397 have a much more 
centrally concentrated spatial distribution than do either 
the fainter CVs or the ABs. Indeed, five of six bright CVs 
are located within about 2 r c of the cluster center, while 
all but one of the fainter CVs lie outside of 4r c . This 
suggests that CVs are preferentially formed by dynamical 
processes within the core and surrounding region and 
migrate to larger radii as they age and undergo repeated 
scattering interactions that act to increase the typical 
size of their orbits in the cluster potential. Since repeated 
strong interactions between compact binaries and other 
stars ultimately result in the ejection of the binaries from 
clusters, the lifetime of a compact binary in a cluster is 
extended if the amount of time it spends in the central 
region of the cluster is reduced. Thus, it is reasonable 
to expect that long-lived CVs in globular clusters should 
generally be found outside of the densest regions. This 
is consistent with the spatial distribution of the CVs in 
NGC 6397. 
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